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ABSTRACT 

Far-infrared and optical [O in] lines are useful tcmeprature-dcnsity diagnostics of nebu- 
lar as well as dust obscured astrophysical sources. Fine structure transitions among the 
ground state levels ls 2 2s 2 2p 3 3 -Po,i,2 give rise to the 52 and 88 fim lines, whereas tran- 
sitions among the 3 Po,i,2, , -D2, 1 Sq levels yield the well-known optical lines AA 4363, 
4959 and 5007 A. These lines are excited primarily by electron impact excitation. But 
despite their importance in nebular diagnostics collision strengths for the associated 
fine structure transitions have not been computed taking full account of relativistic 
effects. We present Breit-Pauli R-matrix calculations for the collision strengths with 
highly resolved resonance structures. We find significant differences of up to 20% in 
the Maxwellian averaged rate coefficients from previous works. We also tabulate these 
to lower temperatures down to 100 K to enable determination of physical conditions in 
cold dusty environments such photo-dissociation regions and ultra-luminous infrared 
galaxies observed with the Herschel space observatory. We also examine the effect of 
improved collision strengths on temperature and density sensitive line ratios. 
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Processes - Atomic Data 



INTRODUCTION 



03 optical lines have long been standard nebular temper- 
ature diagnostics with wavelengths almost in the middle of 
the optical spectrum at AA 4363, 4959, 5007 (viz. Aller 1956, 
Dopita and Sutherland 2003, Pradhan and Nahar 2011). In 
recent years, owing to the advent of far-infrared (FIR) space 
observatories and instruments such as the Intfrared Space 
Observatory - Long Wavelength Spectrograph (ISO-LWS), 
the Spitzer Infrared Spectrograph, and the Herschel Pho- 
todetector Array Camera and Spectrometer (PACS) the 03 
FIR lines have proven to have great potential in providing di- 
agnostics of physical conditions in a variety of astrophysical 
objects that are generally obscured by dust extinction at op- 
tical or shorter wavelengths. These range from Galactic H n 
regions (Martin-Hernandez etal. 2002, Morisset etal. 2002) 
to star- forming galaxies at intermediate redhift (Liu etal. 
2008) and ultra-luminous infrared galaxies (ULIRGs). For 
example, the 03 FIR lines at AA 88 and 52 /im are observed 
from dusty ULIRGs, which are copious IR emitters and be- 
come more prominent with increasing redshift (Nagao et al. 
2011). They may be valuable indicators of the metallicity 
evolution from otherwise inaccessible star-forming regions 
buried deep within the galaxies (Houck etal. 2004, 2005). 

The forbidden FIR lines arise from very low-energy ex- 



citations within the fine structure levels of the ground state 
of atomic ions, such as the 03 3 P — > 3 Pi transition at 88.36 
]iim and the 3 Pi — > 3 P2 transition at 51.81 /an . As such 
they can be excited by electron impact at low temperatures, 
even T e ~ 1000 K or less. That also accounts for their utility 
since the FIR lines can be formed in (and therefore probe) 
not only H 11 regions but also photo-dissociation regions 
(PDRs) where the temperaure-density gradients are large 
(Nagao etal. 2011). 

However, excitation of levels lying very close to each 
other implies that the associated cross sections need to be 
computed with great accuracy at very low energies in or- 
der to yield reliable rate coefficients. The Maxwellian elec- 
tron distribution at low temperatures samples only the near- 
threshold energies above the small excitation energy of the 
fine structure transition. Relativistic fine structure separa- 
tions therefore assume special importance even for low-Z 
atomic ions in determining not only the energy separation 
but also the interaction of the incident electron with the tar- 
get levels. Owing to its prominence in astrophysical spectra, 
a large number of previous studies have been carried out 
on electron impact excitation of O ill (viz. compilation of 
evaluated data by Pradhan and Zhang 2001). Among the 
recent ones, whose collision strengths are employed in astro- 
physical models, are Burke etal. (1989) and Aggarwal and 
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Keenan (1999). But these calculations are basically in LS 
coupling (Burke etal. 1989), or with intermediate coupling 
effects introduced perturbatively via an algebraic transfor- 
mation from the LS to LSJ scheme (Aggarwal and Keenan 
1999). Although the earlier calculations employed the cou- 
pled channel R-matrix method which takes account of the 
extensive resonance structures, the fine structure separa- 
tions are not considered. In this report we take account of 
both the resonances and fine structure in an ab initio man- 
ner. 

Another recent development in relativistic R-matrix 
codes is the inclusion of the two-body fine structure Breit in- 
teraction terms in the Breit-Pauli hamiltonian (Eissner and 
Chen 2012, Nahar etal. 2011). A relativistic calculation of 
collision strengths can therefore be carried out, including 
fine structure explicitly and more accurately than in previ- 
ous works. Relativistic effects are likely to be insignificant 
for optical transitions compared to the FIR transitions since 
the former involve relatively larger energy separations and 
relatvistic corrections are small. Nevertheless, we consider 
all 10 forbidden transitions among the levels dominated by 
the ground configuration of O III. 



2 THEORY AND COMPUTATIONS 

A brief theoretical description of the calculations is given. 
In particular, we describe relatvistic effects and the repre- 
sentation of the (e + ion) system. 



2.1 Relativistic fine structure 

The relativistic Hamiltonian (Rydberg units) in the Breit- 
Pauli R-matrix (BPRM) approximation is given by 
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where the last three terms are relativistic corrections, 
spectively: 



the mass correction term, // mass — _ 
the Darwin term, if Dar = ^f- £ . V 
the spin — orbit interaction term, 
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Eq. (2) representes the one-body terms of the Breit in- 
teraction. In addition, another version of BPRM codes in- 
cluding the two-body terms of the Breit-interaction (Nahar 
etal. 2011; W. Eissner and G. X. Chen, in preparation) has 
been developed, and is employed in the present work. 



2.2 Effective collision strengths 

Cross sections or collision strengths at very low energies may 
be inordinately influenced by near-threshold resonances. 
Those, in turn, affect the effective collision strengths or rate 
coefficients computed by convolving the collision strengths 
over a Maxwellian function at a given temperature T as 



T ij (T e )= / n ij (e)exp(-e/fcT e )d(e/fcT e ), 



(3) 
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Figure 1. Collision strengths for the 03 IR fine structure tran- 
sitions 2p 2 ( 3 P - 3 Pi, 3 Pi - 3 P2) at A A a) 88 (im and b) 52 (im 
respectively. High resolution at near-threshold energies is neces- 
sary for accuracy in rate coefficients at low temperatures. The 
top panel shows an expanded view in the region E C 1 Rydberg; 
both transitions have similar resonance structures. 
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Figure 2. Maxwellian averaged effective collision strengths T(T) 
(Eq. 1) for the transitions q — 3 Pi at 88 /im and 3 Pi — 3 P2 at 52 
fim (solid lines, c.f. Fig. [TJ. Previous results without relativistic 
(l)cffccts (Aggarwal and Keenan 1999) are also shown (dashed lines) 
in the temperature range available T e ^ 2500 K. 



(2) 



where Eij is the energy difference and Qij is the col- 
lision strength for the transition i j. The exponentially 
decaying Maxwellian factor implies that at low temperatures 
only the very low energy Qij(E) would determine the T(T). 



2.3 Wavefunction representation and calculations 

Based on the coupled channel approximaton, the R- 
matrix method (Burke etal. 1971) entails a wavefunc- 
tion expansion of the (e + ion) system in terms of 
the eigenfuctions for the target ion. In the present case 
we are interested in low- lying O/FIR transitions of the 
ground configuration 2s 2 2p 2 . Therefore we confine our- 
selves to an accurate wavefunction representation for the 
first 19 levels dominated by the spectroscopic configurations 
[ls 2 ]2s 2 2p 2 , 2s2p 3 , 2s 2 2p3s. A much larger set of correlation 
configurations is included for configuration interaction with 
the spectroscopic terms using the atomic structure code SU- 
PERSTRUCTURE (Eissner etal. 1974, Nahar etal. 2003): 
[ls 2 ]2p 4 , 2s 2 2p3p, 2s 2 2p3d, 2s 2 2p4s, 2s 2 2p4p, 2s2p 2 3s, 
2s2p 2 3p, 2s2p 2 3d, 2s 2 3s 2 , 2s 2 3p 2 , 2s 2 3d 2 , 2s 2 4s 2 , 2s 2 4p 2 , 
2s 2 3s3p, 2s 2 3s4s, 2s 2 3p3d, 2p 3 3s, 2p 3 3p, 2p 3 3d. We note here 
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that the crucial fine structure separations between the 
ground state 3 Po,i.2 levels reproduced theoretically agree 
with experimentally measured values to ~3% (Nahar etal. 
2011, see Pradhan and Nahar (2011) for a general descrip- 
tion of atomic processes and calculations). The observed 
energies were substituted for theoretical ones in order to 
reproduce the threshold energies more accurately. This is 
of particular importance for excitation at low temperatures 
dominated by near-threshold resonances. Even though the 
observed and experimental values are close, a small differ- 
ence in resonance positions relative to threshold can introuce 
a significant uncertainty in the effective collision strengths. 

The collision strengths were computed employing the 
extended Breit-Pauli R-matrix (BPRM) codes (Eissner and 
Chen 2012). Particular care is taken to test and ensure con- 
vergence of collision strengths with respect to partial waves 
and energy resolution. Total (e + ion) symmetries up to 
(LS)J7T with J ^ 19.5 were included in the calculations, 
though it was found that the collision strengths for all for- 
bidden transition transitions converged for J ^ 9.5. An en- 
ergy mesh of AE < 10~ 4 Rydbergs was used to resolve the 
near-thresold resonances. The resonances were delineated in 
detail prior to averaging over the Maxwellian distribution. 

3 RESULTS AND DISCUSSION 

We describe the two main sets of results for the FIR and the 
optical transitions, as well as diagnostics line ratios. 

3.1 Far-infrared transitions 

The BPRM collision strengths for the two FIR fine struc- 
ture transitions AA 88, 52 /im are shown in Fig. la,b. Al- 
though the resonance structures look similar the magnitude 
and energy variations are not the same. The Maxwellian av- 
eraged effective collision strengths T(T) are quite different, 
as shown in Fig. 2. While T(T; 3 P - 3 Pi) for the 88 ^m tran- 
sition is relatively constant over three orders of magnitude 
in temperature, the T(T; 3 Pi — 3 P2) for the 52 /im transi- 
tion varies by about a factor of 1.5 from the low-temperature 
limit of 100 K to temperatures T > 10,000 K. A compari- 
son with the earlier work by Aggarwal and Keenan (1999) 
is shown as dashed lines, which range down to their lowest 
tabulated temperature 2500K. It can be noted that if the Ag- 
garwal and Keenan values are extrapolated linearly to lower 
temperatures then one would obtain fairly constant effective 
collision strengths. But the present results show marked dif- 
ference owing to resonance structures as in Fig. 1. 

Such temperature sensitivity of the otherwise density 
sensitive 52/88 line ratio is illustrated in Fig. 3a, b. In Fig. 3a 
the solid lines are ratios with the present collision strengths, 
and the dashed lines are using previous results (Aggarwal 
and Keenan 1999). We find very good agreement, imply- 
ing that at all temperatures down to 2500K the differences 
in line ratios would be negligble. Fig. 3b shows the 52/88 
ratio at various temperatures between 100K and 10,000K. 
Whereas the ratio is relatively constant with density at 
100K, its dependence on density varies significantly with 
increasing temperature. The density-temperature dianostic 
value of the 52/88 ratio is apparent from these curves. 
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Figure 3. The density and temperature dependence of the AA 
52/88 fim line ratios; a) the solid curves are present results at 
2500K and 30000 K and the dashed lines are using the Aggar- 
wal Keenan (1999) effective collision strengths; b) 52/88 ratio at 
temperatures 100 < T(K) ^ 10000. 

Therefore, care must be exercised to establish a tempera- 
ture regime for the emitting region. Fig. 3a shows line ratios 
computed at 2500 K and 10,000K, and variations with elec- 
tron density. It is also found that the values of line ratios at 
2500K and 1000K are very close together, implyling cover- 
gence for T ~ 1000K. Fig. 3b clearly demonstrates that the 
line ratio decreases rapidly for T < 1000K to almost flat at 
100K. The low-temperature regime 100-1000K is therefore 
indicated by the curves shown in Fig. 3b, as well as the limit 
where the 52/88 ratio is temperature invariant. So the 52/88 
ratio is excellent density diagnostics in the typical density 
range Log N e ~ 3-4 for T > 1000K without much depen- 
dence on temperature (Fig. 3a). However, at lower temepra- 
tures the ratio may differ by up to a factor of ten (Fig. 
3b). Whereas the primary variations are owing to the expo- 
nential factors in T(T), (Eq. 3), we emphasize the role of 
relativistic fine structure splitting between the 3 Po,i,2 levels 
and near-threshold resonances lying in between. 

3.2 Forbidden optical transitions 

Fig. 4 shows the collision strengths for the optical transi- 
tions 3 Pi - 1 D 2 , 3 P 2 - 1 D2, 1 So - 1 D 2 at A A 4959, 5007, 4363 
respectively. The effective collision strengths for the 03 op- 
tical lines are shown in Fig. 5. These also differ significantly 
from previously available ones, by up to 15%. The new re- 
sults are also obtained down to 100K; their limiting values 
at low temperatures tend to 0.4:0.6:1.0. Since AA 4959, 5007 
are often blended, it is common to plot the bleneded line 
ratio (4959+5007)/4363 shown in Fig. 6. This ratio varies 
over orders of magnitude since the upper-most So level is 
far less excited at low energies than the 1 L>2, and therefore 
the level populations and line intensities depend drastically 
on temperature. A comparison is made with fine structure 
collision strengths derived from the LS term values of Ag- 
garwal and Keenan (1999) divided according to statistical 
weights, again shown as dashed lines in Figs. 5 and 6. How- 
ever, similar to Fig. 3a), the differences in effective collision 
strengths does not translate into any significant differences 
in line ratios even at Log T = 4.5 (« 30000K) wwhere the 
values differ most. 
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Figure 4. Collision strengths of the 03 optical transitions. The 
two transitions 1 D 2 — 3 ft, 1 D 2 — 3 P2 in a) and b) have similar 
resonance structures; the top panel a) presents an expanded view 
below E < 1 Rydberg. 
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Figure 5. Effective collision strengths of the 03 optical transi- 
tions 1 D 2 - 3 Pi, 1 D 2 - 3 P2, 1 So - 1 D 2 at AA 4959, 5007 and 4363 
respectively (c.f. Fig. 4). 
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Figure 6. Blended 03 line ratio (4959 + 5007)/4363 vs. tem- 
perature, at N e = 10 3 cm -3 . The dashed line using earlier data 
(Aggarwal and Keenan 1999) is plotted down to 2500K. 



3.3 Maxwellian averaged collision strengths 

In Table 1 we present the effective collision strengths (Eq. 
3) for the 10 transitions among the ground configuration 
levels and their wavelengths. The tabulation is carried out 
at a range of temperatures typical of nebular environments, 
including the low temperature range T ^ 1000K not hereto- 
fore considered. 



3.4 Conclusion 

Improved collision strengths including fine structure with 
relativistic effects are computed. Owing to the diagnostic 
importance of the 03 forbidden FIR and optical lines, the rel- 
atively small but significant differences of up to 20% should 
provide more accurate line ratios. Particular attention is 
paid to the resolution of resonances in the very small en- 
ergy region above threshold(s), enabling the study of low 
temperature behavior. 

The line emissivities and ratios computed in this work 
demonstrate the temperature-density behaviour at low tem- 
peratures and at typical nebular temperatures. However, 
depending on the astrophysical sources a complete model 
of line emissivities may also need to take into considera- 
tion the Bowen fluorescence mechanism: the radiative ex- 
citation of 2p 2 3 P 2 - 2p3d 3 P 2 ° by He 11 Lya at 304 A and 
cascades into the upper levels of the forbidden transitions 
considered herein (viz. Saraph and Seaton 1980, Pradhan 
and Nahar 2011). In addition, for higher temperatures T > 
20,000K proton impact excitation of the ground state fine 
structure levels 3 Po,i,2 needs to be taken into account; at 
lower temperatures the excitation rate coefficent due to elec- 
trons far exceeds that due to protons (Ryans etal. 1999). 
Finally, there may be some contribution from (e + ion) re- 
combination from O IV to 03, since recombination rate co- 
efficients increase sharply towards lower temperatures while 
collisional excitation rates decrease (level-specific and to- 
tal recombination rate coefficients may be obtained from 
S. N. Nahar's database NORAD at: www.astronomy.ohio- 
state.edu$simna/iar). Recombination contribution depends 
on the O iv/03 ionization fraction, which at low tempera- 
tures would be small. 
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